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This  report  presents  the  results  of  a study  on  pulse  applications 
I and  parameter  identification  for  the  U.5.  Navy  medium-weight  shock  machine. 

G.J.  O'Hara  v<as  project  manager  for  the  Naval  Research  Laboratory. 
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information  in  the  loading  effects  of  equipment  on  structural  motions. 
Professor  S.F.  Masri,  University  of  Southern  California,  performed  the 
investigations  on  parameter  identification. 
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PULSE  APPLICATIONS  AND  PARAMETER  IDENTIFICATION 
FOR  U.S.  NAVY  MEDIUM-WEIGHT  SHOCK  MACHINE 

SECTION  1 
INTRODUCTION 


1 . 1 BACKGROUND 

The  U.S.  Navy  medium-weight  shock  machine  has  played  a paramount 
roie  over  the  past  30  years  in  the  hardening  of  shipboard  equipment  to  shock 
ioads.  Shock  ioading  of  shipboard  equipment  originates  from  depth  charges, 
conventionai  and  nuclear  bombs,  torpedoes,  and  missiles.  Supplementing  the 
shock  machines,  shock  barges  were  subsequently  placed  into  service  by  which 
equipment  mounted  in  the  barges  were  subjected  to  underwater  detonations. 

The  combination  of  the  shock  test  machines  and  the  shock  barge 
has  proved  singularly  effective  in  the  upgrading  and  qualification  of  equipment 
to  specified  hardness  levels  in  a broad  sense.  However,  rapid  changes  in 
naval  warfare  require  more  refined  test/analyses  and  continuous  improvements 
thereof.  These  improvements  are  necessary  to  meet  the  ever-increasing  threats 
caused  by  technical  change,  particularly  those  due  to  the  improving  detection 
and  accuracy  of  offensive  weapons  systems.  Balanced  hardness  in  a navat 
weapons  system  is  becoming  a more  critical  requirement  to  enhance  survivability 
and  operability,  both  during  and  after  an  attack. 

Failures  and  malfunctions  of  both  equipment  and  structure  closely 
correlate  to  threat  time-histories.  This  correlation  motivates  the  continuing 
need  for  test  simulation  ever  more  closely  approximating  threat  environments. 

The  invention  and  development  of  a high-force  output  mechanical  pulse  generator 
(Refs.  1,  2)  showed  considerable  promise  in  inducing  structural  motions  simu- 
lating those  motions  induced  by  threats.  The  use  of  a mechanical-force  pulse 
generator  with  the  U.S.  Navy  medium-weight  shock  machine  should  extend  test 
machine  applications  and  provide  for  a more  flexible  transient  shock  capability. 

Another  critical  area  in  need  of  development  Is  an  accurate 
description  of  the  structural  load  paths  between  input  threat  and  the  mounting 
locations  of  operational  and  weapons-system  equipment.  With  an  adequate 
Not*:  Manuacript  lubmittad  Nov*mb*r  7,  1077, 
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structural  description,  the  motions  at  equipment  mounting  generated  by  an 
attack  may  be  predicted  and  the  structural  path  can  be  analyzed  for  methods 
mitigating  shock  transmission.  Obtaining  the  appropriate  transmitted  shock 
time-histories  at  the  mounting  points  of  equipment  permits  hardening  of  the 
equipment.  One  of  the  better  and  more  practical  descriptions  of  the  structural 
load  paths  is  impedance  and  mobility.  The  severity  of  postulated  attacks 
will  range  from  linear  into  varying  degrees  of  the  nonlinear  regions  such 
that  high-force-level  impedance  measurements  are  required.  The  high-force 
levels  are  required  for  impedance  measurements  in  a quasi-1 inear  sense  or 
for  the  formation  of  nonlinear  functionals.  Electrodynamic  shakers  are 
available  with  high-force  outputs,  and  the  mechanical  force  pulse  generator 
discussed  in  the  previous  paragraph  may  also  be  used. 

Impedance  measurements  may  be  used  to  extract  the  normal  modes, 
damping,  and  resonant  frequencies  of  a structure.  Another  approach  to  handling 
impedance  measurements  is  the  conversion  from  their  nonparametr ic  form  to 
parametric  form.  In  this  latter  procedure,  parameters  are  identified  for  each 
structural  path  measured,  and  represent  a coupled  and  distributed  system. 
Representation  of  impedance  in  parametric  form  provides  additional  opportu- 
nities for  design  changes  or  modifications  for  shock  transmission  mitigation. 


1 .2 


SCOPE 


The  scope  of  this  study  covers  the  following  two  objectives: 


• Provide  the  design  for  a mechanical -force  pulse  generator, 
specify  associated  hydraulic  power  unit  and  controls,  and 
provide  pulse-train  profiles  and  mandrels^  The  force-pulse 
system  will  be  designed  for  use  with  the  U.S.  Navy  medium- 
weight  shock  machines  and  for  use  in  measurements  of  impedance. 


• -^Perform  an  exploratory  parametric  identification  study  of 

U.S.  Navy  medium-weight  shock  machine,  test  articie,  and  impact 
loads  using  NRL-furnished  data. yK  Determine  equivalent  system 
model  and  associated  parameters 
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SECTION  2 


MECHANICAL  FORCE-PULSE  GENERATOR 


2.1  REQUIREMENTS 

Transient  shock  tests  on  equipment  and  systems  to  simulate  the 
motions  induced  by  a conventional  explosive  or  nuclear  attack  are  largely 
limited  to  single-axis  test  machines.  Further  limitations  exist  in  the  size 
and  weigtit  of  equipment  that  can  be  tested.  Simulating  multiaxis  loading 
on  large  equipment  with  many  degrees  of  freedom  represents  a difficult  problem, 
as  it  is  impractical  to  generate  continuously  varying  forces  of  sufficient 
magnitude.  This  problem  becomes  extremely  difficult  or  impossible  where 
in-place  or  field  tests  are  required.  On  the  other  hand,  short  duration 
forces  of  large  magnitudes  over  a wide  frequency  range  can  be  generated  by 
mechanical  pulse  generators.  Since  a discrete  number  of  pulses  superficially 
presents  an  appearance  quite  different  from  a continuous  excitation  signal, 
it  becomes  necessary  to  select  the  pulses  in  such  a way  that  *the  resulting 
vibration  of  the  structure  matches  as  closely  as  possible  the  response  (e.g., 
displacement,  velocity,  or  acceleration)  produced  by  the  continuous  force, 
as  determined  by  an  appropriate  error  criterion.  This  approach  is  shown 
i n Figure  1 . 

It  is  important  to  note  that  the  method  of  Figure  1 requires  that 
the  criterion  response  to  the  continuous  input  be  known,  which  generally  would 
not  be  true  in  practice.  To  accomplish  this  objective,  the  approach  proposed 
here  assumes  that:  (l)  a mathematical  model  or  impedance  measurements  of  the 
system  under  study  are  known  and  (2)  the  inputs  of  interest  (e.g.,  explosive 
or  nuclear  blast)  are  given.  Under  these  conditions  the  "criterion  response" 
can  be  calculated  and  used  to  obtain  the  pulse  train  for  the  simulated  test. 

In  general,  the  response-time  history  of  a test  article  under  simu- 
lated test  should  show  a reasonable  approximation  to  the  expected  environmental 
phenomena  for  meaningful  hardness/vulnerability  evaluation.  To  accomplish 
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Optimum  rccpoiue  simulation  of  multidegrce  systems  by  pulse  excitation 


the  desired  response  motion,  a computer  algorithm  is  used  for  the  optimum 
selection  of  a finite  number  of  pulse  heights,  duration,  and  onset  times  to 
accomplish  this  objective. 

The  basic  criterion  used  is  the  integral  squared  error  between  the 
reference  and  simulated  response,  evaluated  at  a sufficient  number  of  points 
within  the  system  under  test  to  characterize  it  as  completely  as  possible. 
Given  the  error  criterion,  the  pulse  occurrence  times,  pulse  widths,  and  the 
pulse  amplitudes  are  selected  by  a systematic  search  algorithm  such  that  the 
error  is  minimized. 

Application  of  the  pulse  generator  for  impedance  measurements 
requires  a different  form  of  the  pulse  train.  Use  of  the  pulse  generator  in 
this  manner  requires  a sufficiently  long  pulse  train  duration  (T)  to  provide 
the  required  frequency  resolution  (f  ■ y")  ^nd  force  levels  sufficient  to 
provide  favorable  s igna I -to-no i se  ratios  over  the  frequency  range  of  interest. 
Generally,  configuration  of  the  pulse  train  is  optimized  to  generate  an 
approximately  flat  Fourier  magnitude  spectrum.  Studies  have  also  shown  that 
diagnostic  pulse  trains  can  be  generated  that  approximate  (in  a transient 
sense)  sine  wave  bursts,  chirp,  and  random  effects.  Both  chirp  and  random 
are  particularly  useful  in  developing  functionals  for  nonlinear  systems. 

The  pulse  generator  inherently  possesses  several  desirable  features 
for  transient  load  application  with  large  systems  when  tested  in  place  and/or 
in  the  field.  These  features  include: 

• High  force  output 

• Wide  frequency  band  from  static 

• Portabi 1 i ty 


• Ganged  operations 

• Design  simpi ici ty 

• Rel iabi I i ty 


2.2  PRINCIPLE  OF  OPERATION 


The  concept  of  a mechanical  pulse  generator  simply  turns  a device 
for  energy  absorption  around  to  obtain  the  desired  force  profiles  as  a useful 
output  (Ref.  3).  By  drawing  a metal  bar  or  mandrel  through  a cutting  tool 
(or  vice  versa)  with  suitable  motive  power  (air  pressure,  hydraulic  pressure, 
explosive  force,  electric,  mechanical),  a series  or  a set  of  force-time 
histories  may  be  generated.  Reaction  at  the  attach  point  of  the  device  trans- 
mits a force  output  to  a test  article.  Figure  2 illustrates  the  device  with 
motive  power  for  operation  supplied  by  the  stored  energy  in  a oneumatic 
cylinder.  Figure  3 is  an  example  of  a force-pulse  time  history  as  output 
from  the  device.  The  device  may  be  started  or  stopped  by  a variety  of  devices. 

Forces  may  be  generated  singly  and  in  series  as  for  one  cutting 
head  or  in  parallel  with  multiple  cutting  heads.  A single  large  cutting  tool 
may  be  used  if  economic.  Cutting  of  metal  may  be  groove  cuts,  l.e.,  the  work 
is  wider  than  the  tool  width;  or  single  cut,  i.e.,  the  tool  is  wider  than  the 
work.  The  cutting  tool  may  have  a wide  variety  of  shapes  to  suit  the  force 
waveform  output  requirement  needed.  Large  forces  may  be  generated  from  this 
device  and  the  force  required  to  cut  metal  is  largely  independent  of  rate 
(velocity).  The  force  is  a function  of  volume  of  chips  cut  (depth,  width,  and 
length  of  cut)  and  the  specific  energy  of  cutting.  The  energy  absorbed  in 
metal  cutting  is  given  by  (Ref.  3) 

FI  “ Iwtu  (l) 

where 

F “ Force  of  cutting,  lb 

1 ••  Length  of  cut.  In. 

w ” Width  of  cut,  in. 

t ■ Depth  of  cut,  in. 

u ■ Specific  energy  of  cutting,  in-lb/cu  in. 
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FORCE,  LB 


Fig.  2 — Pulse-forming  device  witK  driving  end  control  system 


I . 

The  values  of  specific  energy  of  cutting  for  mild  steel  and  aluminum  are 
3 X 10  in.-lb/cu  in.  and  1.5  x 10^  in.-lb/cu  in.,  respectively.  A load  cell 
or  strain  gage  may  be  incorporated  in  series  with  the  device  to  provide  force- 
time history  readout  as  the  device  is  operated. 

Metal-cutting  information  from  the  Metals  Handbook  (Ref.  A)  generally 
supports  the  foregoing  information.  Higher  rake  angles  (to  20  deg)  on  the 
cutting  tool  tend  to  reduce  cutting  force  and  data  scatter.  At  low  cutting 
speeds  aluminum  is  rate  sensitive,  with  the  cutting  force  changing  somewhat 
exponentially  between  AO  ips  to  100  ips.  It  is  relatively  constant  at  higher 
speeds.  Nominal  power  requirements  for  cutting  various  metals  are  given  as 
fol lows: 


Magnes i urn 

0.10 

hp/in.^/min 

Aluminum 

0.15 

hp/in.^/min 

Copper  a I loys 

0.25 

hp/in.^/min 

Steels 

0.80 

hp/in.  Vmin 

These  values  may  change  considerably  with  hardness  and  alloy 
content.  Calibration  is  required  for  each  mandrel  material  selected. 

2.3  PULSE  DEVICES 

Several  pulse  devices  have  been  constructed  and  successfully 
operated.  The  original  experimental  system  is  shown  in  Figures  ^4  and  5 
(Ref.  1).  The  device  in  Figure  6 was  used  to  simulate  response  motions  in 
large  equipment  induced  by  a nuclear  attack  (Refs.  5.  6,  7). 

Comparison  of  the  predicted  response  (determined  from  impedance 
measurements),  pulse-simulated  response,  and  measured  response  are  provided 
in  Figure  7<  This  same  pulse  generator,  which  was  also  used  to  measure 
transfer  impedances  (Ref.  8),  has  a nominal  output  force  capacity  to 
10,000  lb.  Another  version  in  the  form  of  drop  shock  tower  is  pictured  in 
Figure  8 and  output  data  is  given  in  Figure  9.  Figure  9 data  are  from  a 
re i nforced-concrete  protective  structure,  flush  buried  to  ground  surface 
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Fi){.  5 — Cutting  fixture  with  mandrel  tool  holder,  right;  bearing  guide 
plate  on  mandrel,  center;  and  output  plate,  left 
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— Control  room  platform:  accalaration-tlma  hiatoriaa  of 
predicted,  pulie-aimulated,  and  pulae-taated  motiona 
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INPUT  ACCELERATION  OF 
DROP  CARRIAGE  OF 
PULSE  GENERATOR 


F>8>  9 — Oscillograph  records  amplitude  histories  of  drop  carriage  and  EDC  structure  response 


where  the  generator  was  configured  for  three  5000-lb  pulses  of  25  msec 
durations  each.  This  drop  shock  tower  has  the  potential  capacity  to  produce 
pulse  amplitudes  up  to  150,000  1b. 

2.^*  NRL  PULSE  GENERATING  SYSTEM 

The  NRL  Pulse  Generating  System,  Figure  10,  was  primarily  designed 
for  laboratory  use  in  conjunction  with  the  medium-weight  shock  machine. 

However,  this  system  is  inherently  flexible  and  may  be  used  In  remote  loca- 
tions and  on  shipboard.  A full  technical  and  engineering  description  of  this 
system  was  provided  to  NRL  in  Report  R-7710-3"^333 , "Specification  for 
Evaluation  of  Impedance  Analysis  Techniques  for  Naval  Applications," 

December  1976  by  Agbabian  Associates,  El  Segundo,  Calif.  Major  design 
change  over  previous  versions  was  in  reconfiguration  so  as  to  reduce  the 
weight  below  gage  (load  cell).  Maximum  capacity  of  this  unit  is  10,000  lb 
force.  To  obtain  maximum  capacity,  a dual  cutter  would  be  required,  and 
bottled  compressed  air  would  be  needed  in  place  of  plant  air. 

The  calibration  mandrel  for  the  system  is  drawn  in  Figure  11,  and 
its  Fourier  magnitude  spectrum  is  given  by  Figure  12.  A three-pulse  mandrel 
Fourier  spectrum  is  shown  in  Figure  13  only  to  present  the  wide  spectrum 
that  may  be  obtained. 

The  pulse-generator  system  is  depicted  in  Figure  lA.  This  system 
operates  from  plant  air  to  an  air/hydraulic  multiplier,  thence  to  a hydraulic 
drive  to  operate  the  pulser.  Orifice  flow  control  gives  cutting  velocities 
from  20  in. /sec  to  136  in. /sec. 

Initial  applications  for  this  device  are  planned  for  the  U.S.  Navy 
medium-weight  shock  machine,  illustrated  in  Figure  15  with  its  impact  hammer. 
Single-pulse  impact  of  this  hammer  Is  In  the  order  of  1 msec.  Whereas,  for 
example,  the  mandrel  In  Figure  11  generates  12  pulses  of  0.8  msec  each  for  a 
total  test  record  length  of  76  msec  at  120  in. /sec  cutter  velocity  or  12  puKes 
of  A. 8 msec  each  for  a test  record  of  1*56  msec  at  a 20  in. /sec  cutter  speed. 
This  latter  application  is  schematically  shown  in  Figure  16. 
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(a)  Fourier  magnitude 

Fig.  13  — Frequency  transrorm  of  force-pulse  input  for  3-puUe  train 


Fig.  15  — U.S.  N»vy  high-impact  shock  machine  for  medium-weight  equipment 
with  mounted  teat  structure 
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PULSER 


(a)  Pulse-driven  anvil  (b)  Pulse-driven  anvil  and 

test  structure 

Fi(.  16  — Shock  machine  anvil  with  and  without  teat  structure  using 
pulse-drive  system 


SECTION  3 


U.S.  NAVY  medium-weight  SHOCK  MACHINE 
PARAMETER  IDENTIFICATION 


3. 1 SHOCK  MACHINE  AND  TEST  ARTICLE 

The  vibration  of  a irxjunting  base  or  element  of  structure  is  altered 
when  another  structure  is  mounted  thereto  or  a previously  mounted  structure 
is  replaced  with  a different  component.  Knowledge  of  the  shock  and  vibration 
environments  at  mounting  locations  within  a naval  vessel  is  a necessity  for 
the  determination  of  hardness  and  survivability  of  equipment.  This  is  a 
continuing  problem  due  to  the  continued  changes  in  the  state-of-the-art 
resulting  in  the  upgrade  and  substitution  of  equipment. 


The  environment  at  mounting  points  of  equipment  is  given  by  the 
following  expression  (Ref.  9): 


V'(w) 


where 


V(w) 
V (w) 


V(lu) 
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(2) 


Complex  velocity  of  mounting  point  without  equipment 
Complex  velocity  of  mounting  point  with  equipment  attached 
Driving  point  impedance  of  mounting  point 
Driving  point  impedance  of  equipment 


When  one  equipment  is  to  be  exchanged  for  another  equipment  and 
the  original  velocity  at  the  equipment  mounting  interface  is  known,  the 
response  with  the  substituted  equipment  may  be  predicted  by  (Ref.  9): 


V (w) 


V(w) 


+ Z 


a 


•33 


(3) 


where  the  parameters  of  Equation  3 are  the  same  as  in  Equation  2,  except 
is  the  driving  point  impedance  of  the  substituted  equipment. 
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Other  forms  of  Equations  2 and  3 ("ay  be  used  and  these  basically 
include  the  transfer  impedance  term  from  the  input  load  points  to  the 
equipment  mounting  locations. 

A practical  application  of  the  above  was  conducted  in  I960  (Ref.  10) 
using  the  medium-weight  shock  machine.  Impedance  measurements  were  made 
on  the  anvil  and  upon  a test  structure  as  sketched  in  Figure  17.  The  output 
impedance  of  the  anvil  for  both  magnitude  and  phase  is  in  Figure  13  and  the 
input  impedance  for  the  test  article  is  given  in  Figure  19.  The  anvil  and 
test  article  were  struck  with  the  hammer  as  shown  in  the  conf igurat ion  of 
Figure  15.  The  velocity  response  from  this  test  is  plotted  in  Figure  20  in 
the  form  of  a velocity  shock  spectrum.  A similar  test  was  conducted  upon 
the  anvil  only  and  the  resulting  data  was  modified  by  Equation  2 to  predict 
the  motion  of  the  anvil  if  the  test  article  was  attached  thereto.  This 
prediction  is  also  plotted  in  Figure  20.  While  the  shapes  of  the  two  curves 
are  in  general  conformity,  the  predicted  amplitudes  exceed  actual  up  to  a 
factor  of  nearly  2.  Several  areas  have  been  considered  to  account  for  the 
above  difference  and  the  most  probable  one  Is  the  extremely  low  level  of 
force  used  to  measure  the  impedance  of  the  components.  As  discussed  in 
Section  1 of  this  report,  a sufficient  measurement  force  is  required  in  a 
quasi-1 inear  sense  to  obtain  suitable  values  of  damping  (material,  joints, 
aerodynami c) . 

An  alternative  means  of  measurement  of  a system  is  with  and  without 
load.  Complex  subtraction  of  these  measurements  using  Equation  2 yields 
the  input  impedance  of  the  load.  An  example  of  this  procedure  is  given  in 
Figure  21.  This  example  is  particularly  interesting  due  to  the  extremely 
large  impedance  mismatch  between  the  base  structure  and  the  spring  isolator- 
supported  load  that  was  detected  by  the  impedance  measurements. 

The  foregoing  tests  on  the  shock  machine  are  planned  to  be  repeated 
using  updated  techniques  of  measurement  and  data  processing.  Recent  experi- 
ence indicates  that  very  useful  results  will  be  obtained  by  which  accurate 
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Fig.  18  “ Output  impedance  anvil  (Z22)  medium-weight  shock  machine 
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FREQUENCY  (CPS) 

Pig.  20  — Shock  spectra  velocity  of  anvil  and  teat  article 


(data  record  and  computed) 
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Drtvini;  point  impedance  tests  of  a reinforced'Concrele  structure 
at  attachment  location  of  a shock  isolator 


predictions  can  be  made  on  the  velocity  effects  of  loading  a structure  with 
new  and  changed  equipment.  Once  thCs  procedure  is  established,  application 
to  operating  weapons  systems  may  be  initiated. 

3.2  PARAMETER  IDENTIFICATION 

3.2.1  PROCEDURE 

Impedance  measurements  discussed  previously  are  nonparametric  and 
may  not  often  be  in  a form  most  convenient  to  understand  the  shock  and 
vibration  transmission  path  or  means  to  reduce  transmissions. 

In  the  parametric  identification  approach,  the  mathematical  struc- 
ture of  the  model  is  postulated,  but  its  paratneters  are  not.  Most  of  the 
work  done  so  far  in  structural  system  identification  has  been  considered 
from  the  parameter  estimation  approach.  The  identification  task  in  the 
parametric  model  approach  eventually  reduces  to  a search  in  parameter  space 
where  system  parameters  are  iterated  repeatedly  until  values  are  obtained 
that  meet  a specified  error  criterion.  Among  the  techniques  used  in  this 
approach  are  gradient  or  random  search  methods.  For  a satisfactory  model, 
once  the  parameters  have  been  identified,  methods  can  be  devised  for 
appropriate  insertion  loss  to  mitigate  shock  transmission.  Extending 
parameter  identification  to  nonlinear  systems,  a series  of  impedance  measure- 
ments at  different  force  levels  are  obtained.  This  approach  provides  a group 
of  nonparametric  impedance  curves  from  which  parameters  are  identified. 

From  the  parameters  obtained,  functional  relations  are  set  forth  in  the  model 
as  for  nonlinear  springs  and  damping  to  provide  an  approximate  model  of  the 
nonlinear  transmission  path. 

Random  search  algorithms  for  parameter  optimization  have  been 
widely  applied  and  documented  (Refs.  11  to  17).  They  have  the  advantage  of 
(1)  leading  to  global  solutions  of  nonlinear  systems,  (b)  guaranteed  con- 
vergence, and  (3)  ease  of  computer  Implementation.  On  the  negative  side. 
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random  search  algorithms  may  converge  very  slowly,  particularly  in  criterion 
surfaces  of  high  dimensionality.  The  basic  algorithm  for  minimization  of  a 
criterion  function  J(a),  where  a - (a,,  a,,  a is  a vector  of 

unknown  parameters,  proceeds  as  follows:  (l)  An  initial  parameter  value 
is  estimated  and  J(a°)  is  evaluated;  (2)  trial  points  a*  where 

is  the  given  permissible  region  in  the  m-dimensional  parameter  space,  are 
selected  from  an  appropriate  probability  density  function  defined  over 
(3)  a successful  point  a'*'  is  one  for  which  < J(g').  The 

sequence  J(a')  thus  converges  to  a local  minimum.  Rather  than  using  "pure 
random  search,"  most  algorithms  are  based  on  a "random  creep"  procedure  in 
which  exploratory  steps  are  confined  to  a hypersphere  centered  about  the 
latest  successful  point  a'.  However,  convergence  by  such  procedures  may 
be  extremely  slow,  since  no  allowance  is  made  for  variations  in  the  nature 
of  the  criterion  function  surface  as  the  search  progresses  towards  a minimum. 

Several  procedures  have  been  tried  in  the  past  to  circumvent  slow 
convergence.  The  addition  of  a bias  vector  to  each  iteration  causes  the 
search  to  favor  successful  directions  (Ref.  lA).  Restriction  of  the  search 
to  a directed  hypercone  has  also  been  proposed  (Ref.  16).  However,  such 
schemes  are  highly  directional  and  lose  some  of  the  flexibility  of  the  random 
search  procedure,  as  in  the  case  of  narrow  "canyons"  or  "ridges"  in  the 
criterion  function  surface.  Rastrigin  (Ref.  11)  has  compared  a random 
search  where  each  step  is  random  in  direction  but  fixed  in  length  (fixed- 
step-size  random  search),  with  a fixed-step-size  gradient  search,  and  shov/ed 
that  under  certain  conditions  the  random  search  is  superior.  It  is  intui- 
tively evident,  however,  that  one  would  obtain  even  better  performance  if  the 
step  size  of  the  random  search  procedure  were  optimized  at  each  step  of  the 
iteration.  If  the  steps  are  too  small,  the  average  improvement  per  step  will 
also  be  small  and  convergence  time  will  be  lengthened.  If  the  steps  are 
too  large,  they  may  overshoot  the  minimum;  and  the  probability  of  improvement 
will  again  be  too  small.  Hence,  some  method  of  adapting  the  step  size  to  the 
local  behavior  of  J(oi')  is  indicated. 
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Schumer  and  Steiglitz  (Ref.  15)  have  tested  one  approach  to  adaptive 
step  size  random  search.  They  use  a fixed-step-size  algorithm,  where  the 
trial  vector  Aa'  is  of  length  I and  distributed  uniformly  over  the  hyper- 
sphere of  radius  i whose  center  is  at  the  latest  successful  parameter 
value  a'.  A random  step  of  size  i and  an  incremental  step  of  size  4(1  + a) 
where  0 < a < 1 are  taken  in  the  same  direction.  The  step  size  that  produces 
the  larger  improvement  is  used  as  the  nominal  length  4 for  the  next  iteration 
If  no  improvement  occurs,  4 is  incrementally  reduced.  Thus,  the  algorithm 
is  capable  of  adapting  its  step  size  as  a function  of  the  search. 

The  algorithm  used  for  the  medium-weight  shock  machine  is  another 
approach  to  the  determination  of  the  optimal  step  size  (Ref.  18).  Rather 
than  a fixed-length  step,  steps  used  are  random  in  both  length  and  direction. 
Hence,  the  adaptation  described  below  is  based  on  the  selection  of  the 
optimal  variance  of  the  step-size  distribution  as  the  search  progresses. 

Large  variances  are  desirable  in  the  early,  exploratory  portions  of  a search. 
However,  in  the  vicinity  of  a local  optimum,  a smaller  value  of  the  standard 
deviation,  a,  will  decrease  the  probability  of  overshoot. 

The  algorithm  for  the  adaptive  random  search  consists  of  alternating 
sequences  of  a gloual  random  search  with  a fixed  value  for  the  step  size 
variance  followed  by  searches  for  the  locally  optimal  o. 

figure  22  illustrates  the  adaptive  algorithm  whereby  a very  wide- 
range  search  selects  the  standard  deviation  best  step  size  (a)  for  the 
coarseness  of  the  increments  used,  followed  by  a sequential  precision  search 
of  finer  increments.  As  the  rate  of  convergence  decreases,  a new  precision 
search  is  made,  but  directed  towards  a smaller  step  size.  At  selected 
iteration  intervals,  the  wide-range  search  is  reintroduced  to  prevent 
convergence  to  local  minima. 

The  adaptive  random  search  method  was  used  to  determine  a pulse 
train  (composed  of  8 pulses)  whose  response  spectrum  matches  a criterion 
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! o-ii'i'nbc  spoctrum  as  closely  as  possible  over  a given  frequency  range 

lo  I.  . Since  each  i.idividual  pulse  in  the  train  Is  characterized 

I!!  in  max 

■ 3 independent  parameters  (amplitude,  aj*,  duration,  dj',  initiation  time,  t. 
a total  of  2^1  parameters  are  needed  to  define  the  8 pulses  of  the  desired 
train.  Ncte  that  the  response  spectrum  of  a pulse  train  is  a nonlinear 

function  of  the  pulse  parameters. 

Figure  23  shows  that  superior  results  are  obtained  by  using  the 

adaptive  features  of  the  random  search  method.  While  both  cases  (a)  and  (b) 

started  wi  tin  the  same  standard  deviation,  “ O-l*  ' the  case  (a)  the 

search  was  conducted  for  2000  steps  with  fixeu  o.  When  the  search  for  the 

fnst  standard  deviation,  a is  conducted  at  a rate  of  iteration  frequency  * 

"00,  the  values  of  a,  were  found  to  change  from  0.1  initially  to  0.01 

best 

ar  itsrotion  = .500,  and  0.001  at  iteration  ==  1000. 

it  is  worth  noting  that  the  search  for  the  optimum  solution  of  thii 
o'-oblem  was  performed  without  benefit  of  knowing  the  answer  in  advance. 


3.2.2  PRELIMINARY  STUDIES  ■■  — 

A preliminary  study  v/as  conducted  on  the  shock  machine  anvil  and 
test  article  (Fig.  15)  given  the  objective  function  (shock  spectrum)  of 
Figure  20.  Major  difficulty  occurred  in  a poor  description  of  the  pulse 
inptit  from  the  hammer.  A 1 msec  rectangular  pulse  was  finally  specified 
and  held  fixed  during  all  subsequent  computations.  Shock  spectrum  object! ve 
functions  also  pose  costly  computer  operations  as  each  parameter  iteration 
invoives  conversion  of  trial  model  time-histories  into  shock  spectra  until 
convergence  is  obtained. 

The  results  of  the  preliminary  parametric  study  showed  (for  the 
four-degree-of“f reedom  model  chosen  to  represent  the  distributed  system) 
reasonable  results.  Parameters  for  the  model  are  listed  in  Table  1 and  the 
shock  response  spoctrum  Is  given  in  Figure  2^4.  Results  of  the  spectrum  may 

be  compared  to  Figure  20. 
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TABLE  1.  MODEL  PARAMETERS  OF  ANVIL  AND  TEST  ARTICLE 
(input  impulse  of  1 msec  duration) 


Mass  Mass 

No-  (slugs) 


Damping, 

1 b-sec  Spring  Rate, 

ft  1 b/ft 


1 

2 

3 

k 


36. 

1550. 

18.2 

X 

10 

3.39 

3.31 

10.9 

X 

lo' 

8.78 

79. 

35.9 

X 

lo' 

3-96 

128. 

6.9 

X 

lo' 

In  view  of  the  above,  the  better  procedure  is  to  use  the  impedance  functions 
in  the  form  of  impulse  functions  as- object i ves . This  approach  eliminates 
the  large  uncertainties  of  the  input  pulse  and  reduces  costs  by  eliminating 
shock  spectra  computations. 


36 


DP^PiNG  r 0-00 


100.00  ’50. 00  200.00  250.00  300. 00  350.00  400.00 

FREQUENCY  (HZ) 

:4  — Shock  spectrum  response  of  anvil  and  test  article 
using  identified  parameters 
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ATTN  TECHNICAL  LIBRARY 


COMMANDER 

U S ARMY  MAT  ♦ MECHANICS  RSCH  CTR 
WATERTOWN*  MA  02172 
ATTN  R SHEA 


COMMANDER 

U S ARMY  nuclear  AGENCY 
FORT  BLISS*  TX  79916 
ATTN  TECH  LIB 


COMMANDER 

U S ARMY  WEAPONS  COMMAND 
ROCK  ISLAND  ARSENAL 
ROCK  ISLAND*  IL  61201 

ATTN  TECHNICAL  LIBRARY 


CHIEF  OF  NAVAL  MATERIAL 
navy  DEPT 

WASHINGTON  D.C.  20360 
ATTN  MAT  0323 
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CHIEF  OF  NAVAL  OPERATIONS 
NAVY  DEPARTMENT 
WASHINGTON  D.C.  20330 
ATTN  OP  03FG 
ATTN  OP  985F 


CHIEF  OF  NAVAL  RESEARCH 
DEPARTMENT  OF  THE  NAVY 
ARLINGTON  VA  22217 

ATTN  N PERRONE  CD  474 
ATTN  TECHNICAL  LIBRARY 


OFFICER  IN  CHATTGE 
CIVIL  ENGR  LAB 
NAVAL  CONST.  BATTALION  CTR 
PORT  HUENEME  CA  93041 
ATTN  R OOELLO 
ATTN  TECH  LIB 


COMMANDER 

NAVAL  ELECTRONIC  SYSTEMS  COMMAND 
NAVAL  ELECTRONIC  SYSTEMS  COMMAND  HOS 
WASHINGTON  D.C.  20360 
ATTN  PME  117-21A 


COMMANDER 

NAVAL  FACILITIES  ENGINEERING  COMMAND 

HEADQUARTERS 

WASHINGTON  D.C.  20390 

ATTN  TECHNICAL  LIBRARY 


SUPERINTENDENT 

NAVAL  POSTGRADUATE  SCHOOL 

MONTEREY  CA  93940 

ATTN  CODE  2124  TECH  RPTS  LIBRARIAN 
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DIRECTOR 

NAVAL  RESEARCH  LABRATORY 
WASHINGTON  O.C.  20375 

ATTN  CODE  2027  TECHNICAL  LiB  (3  COPIES) 

ATTN  CODE  8A0  J GREGORY 
ATTN  CODE  84A0  F ROSENTHAL 

ATTN  CODE  8A03A  G OHARA  (4  COPIES) 

ATTN  CODE  8A42  H HUANG 


COMMANDER 

NAVAL  SEA  SYSTEMS  COMMAND 
NAVY  department 
WASHINGTON  0*C*  20362 

ATTN  ORD  - 91313  LIB 
ATTN  CODE  03511  C POHLER 


COMMANDER 

NAVAL  SHIP  RSCH  AND  DEVELOPMENT  CTR 
UNDERWATER  EXPLOSIONS  RSCH  DIVISION 
PORTSMOUTH  V.A.  23709 
ATTN  E PALMER 
ATTN  technical  LIBRARY 


commander 

NAVAL  SHIP  ENGINEERING  CENTER 
CENTER  BUILDING 
HYATTSVILLE  MD  20782 


ATTN 

NSEC 

61200 

ATTN 

NSEC 

6110.01 

ATTN 

NSEC 

6105G 

ATTN 

NSEC 

6105 

ATTN 

6105C1 

ATTN 

TECHNICAL  LIBRARY 

commander 

NAVAL  SHIP  RESEARCH  AND  DEV  CENTER 
BETHESDA  MD  2003A 


ATTN  CODE  17  WW  MURRAY 

ATTN  CODE  142-3  LIBRARY 

ATTN  CODE  174 

ATTN  CODE  11 

ATTN  CODE  2740  Y WANG 

ATTN  CODE  1962 

ATTN  CODE  1903 

ATTN  CODE  1731C 

ATTN  CODE  1171 

ATTN  CODE  19 


commander 

NAVAL  SURFACE  WEAPONS  CENTER 

white  oak 

SILVER  SPRING  MD  20910 

ATTN  CODE  241  J PETES 
ATTN  CODE  1224  NAVY  NUC 
ATTN  CODE  730  TECH  LIB 
ATTN  CODE  240  H SNAY 
ATTN  CODE  243  G YOUNG 


COMMANDER 

NAVAL  SURFACE  WEAPONS  CENTER 
DAHLGREN  LABORATORY 
DAHLGREN  VA  22448 

ATTN  TECHNICAL  LIBRARY 


COMMANDER 

NAVAL  UNDERSEA  CENTER 
SAN  DIEGO*  CA  92152 

ATTN  TECHNICAL  LIBRARY 


COMMANDER 

NAVAL  WEAPONS  CENTER 
CHINA  LAKE  CA  93555 

ATTN  CODE  533  TECH  LIB 


I 


COMMANDING  OFFICER 
NAVAL  WEAPONS  EVALUATION  FACILITY 
KIRTLANO  AIR  FORCE  BASE 
ALBUQUERQUE  NM  87117 

attn  technical  library 


DIRECTOR 

STRATEGIC  SYSTEMS  PROJECTS  OFFICE 
NAVY  DEPARTMENT 
WASHINGTON  0 C 20376 
ATTN  NSP-272 
ATTN  NSP-A3  TECH  LIB 


AF  CAMBRIDGE  RSCH  LABS*  AFSC 
L.G.  HANSCOM  FIELD 
BEDFORD  MA  01730 

ATTN  5UOL  AFCRL  RSCH  LIB 


HEADQUARTERS 

AIR  FORCE  SYSTEMS  COMMAND 
ANDREWS  AFB 
WASHINGTON  0 C 20331 

ATTN  TECHNICAL  LIBRARY 


COMMANDER 

ARMAMENT  DEVELOPMENT6TEST  CENTER 
ELGIN  AFB  FL  325a2 

ATTN  TECHNICAL  LIBRARY 


LOS  ALAMOS  SCIENTIFIC  LAB 

P 0 BOX  1663 

LOS  ALAMOS  NM  8754A 

ATTN  DOC  CONTROL  FOR  REPORTS  LIB 
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SANOIA  LABS 
LIVERMORE  LAB 
P O BOX  969 
LIVERMORE  CA  9A550 

ATTN  DOC  CON  FOR  TECH  LIB 


SANDIA  LABORATORIES 
P.O.  BOX  5800 
ALBUQUERQUE  NM  87115 

ATTN  DOC  CON  FOR  TilAl  SANDIA  RPT  COLL 


U S energy  RSCH  6 DEV  ADMIN 
DIVISION  OF  HEADQUARTERS  SERVICES 
LIBRARY  BRANCH  G-0A3 
WASHINGTON  D C 20545 

ATTN  DOC  CONTROL  FOR  CLASS  TECH  LIB 


UN  IV  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LAB 
P.O*  BOX  808 
LIVERMORE  CA  94550 

ATTN  TECHNICAL  LIBRARY 


AGBABIAN  ASSOCIATES 
250  NORTH  NASH  STREET 
EL  SUGONDO  CA  90245 
ATTN  M AGBABIAN 


BATTELLE  MEMORIAL  INSTITUTE 
505  XING  AVENUE 
COLUMBUS  OH  43201 

ATTN  technical  LIBRARY 
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BELL  TELEPHONE  LABORATORIES  INC. 

MOUNTAIN  AVE 

MURRAY  HILL  NJ  0797A 

ATTN  TECH  RPT  CTR 


BOEING  COMPANY 
P.O.  BOX  3707 
SEATTLE  WA  9812A 

ATTN  AEROSPACE  LAB 


CAMBRIDGE  acoustical  ASSOC 
1033  MASSACHUSETTS  AVE 
CAMBRIDGE  MA  02138 
ATTN  M JUNGER 


( 

( 

CIVIL  / NUCLEAR  SYSTEMS  CORP 
1200  UNIVERSITY  N.  F. 
ALBUQUERQUE  NM  87102 
ATTN  T DUFFY 


electric  boat  DIV 

GENERAL  DYNAMICS  CORP. 
GROTON  CN  063A0 

ATTN  L.  CHEN 


GENERAL  ELECTRIC 
TEMPO-CENTER  FOR 
816  STATE  STREET 
SANTA  BARBARA  CA 
ATTN  DASIAC 


CO. 

ADVANCED  STUDIES 
(P.O.  DRAWER  00) 
93102 


ITT  research  INST 
10  WEST  35TH  ST 
CHICAGO  IL  60616 

ATTN  TECHNICAL  LIBRARY 
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•>i 


INST  FOR  DEFENSE  ANALYSIS 
ADO  ARMY  NAVY  DRIVE 
ARLINGTON  VA  22202 

ATTN  IDA  LIBRARIAN  R SMITH 


J.L.  MERRITT 

CONSULTING  + SPECIAL  ENGR  SVS  INC 
P.O.  BOX  1206 
REDLANDS  CA  92373 

ATTN  TECHNICAL  LIBRARY 


XAMAN  AVIDYNE 

DIV  OF  KAMAN  SCEINCES  CORP 
83  SECOND  AVE 
NW  INDUSTRIAL  PARK 
BURLINGTON  MA  01803 

ATTN  E CRISCIONE 
ATTN  TECHNICAL  LIBRARY 
ATTN  G ZARTARIAN 


KAMAN  SCIENCES  CORP. 

P.O.  BOX  7463 

COLORADO  SPRINGS  CO  80933 

ATTN  technical  LIBRARY 


LOCKHEED  MISSILES  AND  SPACE  CO. 
3251  HANOVER  ST 
PALO  ALTO  CA  94304 

ATTN  TECH  INFO  CTR  D/COLL 
ATTN  T GEERS  D/52-33  BLDG  205 


NATHAN  M.  NEWMARK 
CONSULTING  ENGINEERING  SERVICES 
1114  CIVIL  ENGINEERING  BLDG 
UR6ANA  IL  61801 

ATTN  N NEWMARK 
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POLYTECHNIC  INST  OF  NEW  YORK 
DEPT  OF  applied  MECH 
333  JAY  street 
BROOKLYN  NY  11201 

ATTN  J KLOSNER 


R+D  associates 

P.O.  BOX  3580 

SANTA  MONICA  CA  90A03 

ATTN  TECHNICAL  LIBRARY 


STANFORD  RESEARCH  INST 
333  RAVENSWOOD  AVE 
MENLO  PARK  CA  9A025 

ATTN  SRT  LIB  ROOM  G021 
ATTN  B GASTEN 
ATTN  G A6RAHAMSON 


TETRA  tech  INC. 

630  N ROSEMEAD  BLVD 
PASEDENA  CA  91107 

ATTN  LI-SAN  HWANG 
ATTN  TECH  LIB 


THE  BDM  CORP 
1920  ALINE  AVE 
VIENNA  VA  22180 

ATTN  TECH  lib 


UNlV  OF  MARYLAND 
DEPT  OF  CIVIL  ENGR 
COLLEGE  PARK  MD  207A2 
ATTN  B BERGER 


n 


[ 

I 


URS  RESEARCH  CO, 
155  BOVET  RD. 

SAN  MATEO  CA  9^*402 
ATTN  TECH  LIB 


WEIDLINGER  ASSOCIATES 
3000  SAND  HILL  ROAD 
BUILDING  4 SUITE  245 
MENLO  PARK  CA  94025 

ATTN  J.  ISENBERG 


WEIDLINGER  ASSOCIATES 
HO  EAST  59TH  STREET 
NEW  YORK*  NY  10022 

ATTN  OR.  M.  BARON 


ASSISTANT  CHIEF  FOR  TECHNOLOGY 
OFFICE  OF  NAVAL  RESEARCH 
ARLINGTON*  VIRGINIA  22217 
ATTN  CODE  200 


teledyne  BROWN  ENG. 

MAIL  STOP  44 
300  SPARKMAN  DRIVE 
RESEARCH  PARK 
HUNTSVILLE*  ALABAMA  35807 
ATTN  DR*  MANU  PATEL 


DIRECTOR 

U S ARMY  WATERWAYS  EXPERIMENT  STN 
P.O.  BOX  631 
VICKSBURG  MS  39180 
ATTN  J STRANGE 
ATTN  W FLATHAU 
ATTN  TECH  LIB  (UNCL  ONLY) 
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AF  INSTITUTE  OF  TECHNOLOGY*  AU 
WRIGHT  PATTERSON  AFB»  OH  <*5433 

ATTN  LIB  AFIT  BLDG  640  AREA  BCJNCL  ONLY) 


NKF  ENGINEERING  ASSOCIATES,  INC. 
87  20  GEORGIA  AVENUE 
SILVER  SPRING.  MD.  20910 

ATTN  DR.  R.O.  BELSMEIM 

M&T  CO. 

2130  ARCH  STREET 
PHILADELPHIA,  PA.  19103 

ATTN  B.  McNAIGHT 
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